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Abstract 



The cable equation is a second order, parabolic, partial differential equation 
that describes the evolution of voltage in the dendrite of a neuron. Here 
we look at the various ways in which lambda ( space const ant /variable space 
constant/ length parameter ) is calculated in the cases of linear, passive trans- 
mission as well as nonlinear active transmission. Changes in morphology are 
taken into account by including the case of tapering dendrite, branched den- 
drites, branched dendrites with taper or flare. The case of variable membrane 
resistance and the relationship between input resistance and space constant 
is explored. Finally the reaction diffusion equation governing the diffusion of 
calcium in dendrites and space constant associated with that is described. 
Keywords: cable equation; reaction diffusion equation; lambda ; space con- 
stant; variable space constant; length parameter ; 

1 Introduction 

When a cell is isopotential, the membrane potential is uniform 
at all points of the cell, depending only on time. This is apt 
for describing signalling in a cell body which can be assumed 
to be a sphere [8,21]. However if one needs to look at electrical 
or diffusional properties of dendrites or axons, one needs to ap- 
proximate these by a cylinder ( or more if branched) . Here the 
geometry of the cell plays a role. In this context, space plays 
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a part along with time. Thus a factor called space constant is 
defined which is related to the diameter d of the cell and thus its 
geometry. In this paper, the notion of space constant is explored 
under various conditions : from passive to active, unbranched to 
branched and nontapered to tapered processes. Space constant 
can also be variable due to changes in membrane resistance R m . 
In the case of a passive infinite and semi infinite cable, space 
constant ( A) is the distance where voltage reaches 37% )of 
its original. When electrotonic length L = 1, the voltage will 
decay to only 65% of its original value at x = A. This is due 
to the boundary conditions at x = A. In the infinite cable an 
infinite cylinder with its associated conductance is attached to 
the artificial boundary at x = A and current flows into this 
cylinder, whereas in a finite cylinder with L = 1 the cylinder 
ends at x = A. For a sealed end, there is no current flow and 
conductance at this boundary. Voltage decay is also affected by 
taper or flare in the main cable or its branches [4]. When A is 
large, the spatial decay of input with distance is small and vice 
versa for small A. On the other hand in the case of active or 
sinusoidal stimulation, A w has the dimension of interelectrode 
distance when the peak value of sinusoidal transmembrane dis- 
tance is attenuated to -. 

e 

It is seen that A for both semi infinite and finite cable is directly 
proportional to input resistance. Furthermore A in an active 
cell or with sinusoidal current injection is proportional to \/~Z ( 
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impedance) and (frequency). 
2 Unbranched dendrite 

Linear cable theory assumes that the cable has a uniform 
diameter and leakage resistance. So it can be used in the 
passive case where there is no time or voltage dependent 
conductances. When R mi the specific membrane resistance is 
much larger than i^,the intracellular resistivity, it can be 
safely assumed that all current flow will be one-dimensional 
along the length of the cable. Under such conditions there are 
equations describing the change in voltage across space and 
time along the cable. 



The cable equation is : 

i d 2 v m (x,t) v m (x,t)-v rest , ov r 



dx 2 r m dt 



~\~ C-m Q ^ Iinj{%it) (1) 



where r a is the intracellular resistance per unit length of cable 
with dimensions of ohms per centimeter and r m is the 
membrane resistance of a unit length of fiber measured in units 
of ohms-centimeter. Rewriting this, 



r m d 2 V m {x,t) _ _ OVrn _ j ( A (0 s 

r a ox z ot 
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This leads to , 



2 d 2 V m (x,t) dV m . . , . , . 

A ^ = + *) ~ ^ es * ~ r mIinj{X, t) (3) 



2.1 A DC : 



Adc = (-)^ (4) 



r a 



As r m = and r a = gi 



A OC = (-^ (5) 



where <i is the diameter of the cable or dendrite as the dendrite 
is modeled as a cable. 

= Voexp(^) (6) 

where Vb is a constant or very slowly varying voltage applied to 
one end of a long cable. See Figure 1. 
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2.2 A 



DC taper' 1 



In a tapering dendrite d changes and thus the Xoctaper should 
reflect that. From [14, equation 3 ] 



^DCtaper = ^DC 



r 



1/2 



1 + 4 > 2 



-1/4 



(7) 



where r is the radius of the cylinder at Jj . Jj IS the distance from 
the soma, ro is the radius of the cylinder at x = 0. This can 
be derived as given below. For a dendrite with taper the cable 
equation is [5 equation 3] : 



dV R m r ,d 2 V dV2dr, / ds, , 

V + t — = — — ( 1 )( — ) 

dt 2Rj dx 2 dx r dx dx 



(8) 



ds .dr. l.i . . 

^ = ( 1 + W 5)5 (9) 

Substituting this in equation 8 and multiplying and dividing by 
r gives 



v | dV = R m r r d 2 V dV2dr , ( dr 2 (-1 

dt~ 2R, r { dx 2 dxrdx ){ W } 2 } { } 
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This can be written as 



V + r 



dt 




) (ii) 



From here it can be deduced that : 



(V) 2 = (Anc)^(l + © 2 )^) 



(12) 



Thus 




(13) 



If the tapering dendrite is sufficiently long to be considered semi 
infinite (0 < x < oo) or infinite (— oo < x < oo), the solutions 
can be obtained in closed form. The boundary conditions for 
the semiinfinite case can be any of the following : voltage clamp, 
sealed end, killed end, current injection at one end, natural ter- 
mination, lumped-soma termination. For infinitely large x , V 
remains bounded. Thus lim\ x \.+ 00 \V(x,t)\ < oo,t > 0[28]. 
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2.3 ^DCvariable'' 



Under certain conditions, R m is not uniform but variable across 
the dendrite [11]. Here the cable equation becomes : 



where 



« d 2 V dV 
(\)\x)-^-T m (x)—-V = (14) 



^dc( X ) = xI R ^ (15) 

is the variable space constant. When R m (x) = constant, Xdc(x) = 
Adc 

The electrotonic length of a dendrite is a measure of its length 
in nondimensional terms. It is written as L = j where / is the 
length. 

London et al [11] show that if the total membrane conductance ( 
G m = 1/Rm) is kept fixed for a given cable,then any nonuni- 
formity of G m over the cylinder reduces the electrotonic length 
L. If it is a sealed end cable, any monotonic increase in G m im- 
proves the voltage transfer from input location to the soma. It 
is further improved if it is synaptic input at a distal location as 
opposed to a current source. 
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The input resistance of an infinite cable is the following ratio in 
the limit that the distance between the current passing electrode 
and recording electrode shrinks to zero [9 equations2.14, 2.15] 



= m = n^l (16) 

I%\x) Iq 



This can be written as 



r_rn = r_aX = 

2A 2 2 y J 



For a semi- infinite cable, Roq is double that of an infinite cable 
as there are two semi-infinite cables to one infinite cable [9, 
equation 2.16] 

fl=. = (r B .r„)i=r^ = !f=«^^ (18) 

In multidimensional systems like dendritic trees, cardiac, smooth 
and skeletal muscle, current spread is in more than one dimen- 
sion . Here an effective 'input resistance' is defined as a function 
of electrode separation r. In the steady state for a 2-D case [8 
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equations 5.17,5.18 ] 



For small values of r, 



7 2tt6 r 



For a three dimensional system [8 chapter 5 ] 

v - = A e ^ = + (21) 

7 47rr 47rr v ^R^ ' v 7 

As the dimensionality of the system increases, the input resis- 
tance becomes more and more insensitive to changes in R m . 
In an infinite cable Ri n oc y/R m , for a two dimensional sheet 
Rin oc log(Rm). In a three dimensional structure such as muscle 
tissue R in oc ^^cjj^)- 

2.4 \ AC : 

As R m and C m change with change in voltage in an active cable, 
A also changes. Thus, to get the actual value of A under these 
circumstances, one needs to look at the problem differently.The 
linear cable equation can be solved for measurements with sinu- 
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soidal currents. Here R m has to be replaced by Z m (membrane 
impedance). The input impedance and length become frequency 
dependent. 

The potential distribution in a linear cable when an a.c current 
is injected using intracellular electrodes can be written as [2]: 



V(x,t) = i(J (*)(Z m rO*)ea;p(-^) (22) 



where Ti = r a 

In a.c analysis, Z m is replaced by a complex number so Xac 
cannot have the same physical meaning as Xoc which is a real 
number. So a new A w is defined. 
Along a neurite, 

^ = e-£ (23) 

v o 

where 



^ = \— (24) 



where Z m is the membrane impedance of a unit length of neurite 
and Ti is the resistance of unit length of cytoplasm. If / is large, 
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transmembrane current is entirely capacitative. 



TTL 



(25) 



juC 

where C is the capacitance of unit length of membrane 

Yl ~ e ^+j)VW (26) 



Xuj ~ , ~ y — ~ — (27) 



2 

Substituting, r a = ^ and C = 2naC m where a is the radius 
of the neurite , Ri is cytoplasmic resistivity in Qcm and C m is 
specific membrane capacitance in \iFaradj 'cm 2 in equation 27 



2na 2 I a / a 1 / d 



u27raC m Ri y UJ C m R l y 27r fC m Ri 2 y 7rfR t C m 

(28) 

where d is the diameter of the dendrite. 
Using this equation 22 can be rewritten as : 



V pp = hl pp (Z m r^)exp(-^) (29) 

Z A w 
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where I pp is the peak to peak current at the point of current 
injection. 

Whenever fast changing membrane potentials are encountered 
either as action potentials or as injected pulses, X w will deter- 
mine the electrotonic spread due to the high frequency compo- 
nents of the fast changing voltages. 

Pettersen and Einevoll [13] express \ac(w) for dendritic sticks 
of a finite length. 



f*z\i m (z)\dz 
Jo \i m {z)\dz 



where lu = 2nf is the angular frequency, \i m \ is the amplitude 
of the sinusoidally oscillating current at position z when a sinu- 
soidal current is injected at the soma. Equation 30 reduces to 
the following for a dendrite of infinite length : 

>&(») = ' . 2 J (31) 

where r is the membrane time constant. Xac is thus dependent 
on frequency. It decreases with increasing frequency. Pettersen 
and Einevoll[12] discuss in great length a low pass filtering effect 
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of extracellular potentials as a result of this dependence of X w 
on frequency. 

It is seen that at high frequencies, transmembrane current is 



gation of signals > 5/ TO . Here is given by equation 28. This 
shows that A w oc ^j. This implies that at higher frequencies, 
the Xoj is lower. 

Using another approach, Lindsay and Rosenberg [10] show that 
an active neuron will have space and time constants which reflect 
dynamic biophysical properties of the membrane. The space 
constant A and time constant T give the spatial and temporal 
rates of exponential decay of the total membrane current during 
an action potential. To calculate this one needs to calculate con- 
duction speed of the action potentials which can be done easily 
using extracellular methods. Thus one can get [10, equation 16] 



where 6 is the speed of movement of the action potential train in 
cm /msec, P is the perimeter of the dendrite, A is the crosssec- 
tional area, C m is the capacitance, g a is the axial conductance 
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purely capacitative. At / = |7tt to , R m has no effect on propa- 




(32) 



in mS/cm and g m is the membrane potential during an action 
potential. [10, equation 17] gives : 



A= t /g (33) 



When g m = g m , A = X DC 



2.5 ^ AC taper'- 
2.5.1 linear taper: 

/ x . r/-r 

r[x) = px + tq, p = — - — 

where r(x) is the radius at any given x along the dendrite, ro is 
the radius of dendrite at point and p is the linear taper. 



1. 2r i 
2 nfKiLm 



(2) 2 (X A Ctapered) 2 = 2 ( 7 p ^ ) ( 35 ) 



XACtapered = t( ° + (36) 
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^ACtapered = ^AcO- + ~~j~) 2 (37) 

do 



2.5.2 exponential taper: 

ln(r /rl) 

r[x) = roexp{—p.x),p = 

6 

where r(x) is the radius at any given x along the dendrite, rg 
is the radius of the dendrite at point and p is the exponential 
taper. 

^ACtapered = ~( T^tt)* ( 38 ) 

2 7rJKiC m 



(Z) 2 (\ A Ctapered? = 2 ^ (egp(-p.:c))* (39) 



X A Ctapered = U , ) * (eXff(-p.ff)) 5 (40) 

2 7TjHiC m 



^ACtapered = ^Ac(exp(-p.x))z (41) 
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3 Branched dendrites: 

3.1 Xdc Non tapering : 

This problem can be solved with or without the use of equivalent 
cylinders. 

3.1.1 Without equivalent cylinders: 

Calculate Xpc for each section separately using the equation 5. 

3.1.2 With equivalent cylinders: 

Rail [22,24] showed that dendritic trees could be collapsed into 
a single equivalent cylinder provided they meet the following re- 
quirements [9]: 

1. R m and Ri values are the same in all branches. 

2. All terminals have the same boundary condition. 

3. All terminal branches end at the same electrotonic distance 
L from the origin in the main branch, where L is the sum of 
the Li values from the origin to the distal end of every terminal. 
L corresponds to the total electrotonic length of the equivalent 
cylinder. 

4. At every branch point, infinite input resistances must be 
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matched. If all cables possess the same membrane resistance 
and intracellular resistivity, this implies 

dT = dT + df (42) 

where d Q is the diameter of the parent branch and d\ and d<i are 
the diameters of the daughter branches. This last condition is 
called the d 3 ^ 2 law. 

If these four conditions are met the equivalent cylinder can be 
considered to be a perfect representation of the entire tree pro- 
vided current is injected in the initial terminal. In case there is 
input to any of the daughter branches, an additional constraint 
should be obeyed : 

5. Identical synaptic inputs, whether current injection or con- 
ductance change, must be delivered to all corresponding den- 
dritic locations. 

In such a situation X^c can be calculated from equation 5 where 
the values inserted should be that for the equivalent cylinder. 
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The diameter D of the equivalent cylinder is given by 



D = 



n 



E< /2 



L .?=i 



2/3 



(43) 



D = doi if there is only one dendrite emanating from the soma. 



A = 



4R, 



1/2 



(44) 



3.2 Xdc - Tapering: 

This problem too can be solved with or without the use of equiv- 
alent cylinders. 



3.2.1 Without equivalent cylinders: 

Xdc is calculated for each individual branch using equation 7. 
Here the values of r,r and ^- will varv from branch to branch. 
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3.2.2 With equivalent cylinders: 

Rail [2 2] has shown that under certain conditions, a tapering tree 
can be reduced to a one dimensional cylinder. The condition is 
[8, equation 7.43]: 

1/4 



nr 3/2 



1 + 



dr 
dx 



= constant (45) 



where n is the number of dendritic branches and r the radius of 
all the branch segments. These are functions of actual distance 
x from the soma. The following condition is to be met here 
[8, equation 7.45] 

dA dZ (46) 

where A is the surface area of the dendrites. Z is the electrotonic 
distance. For the dendritic tree equivalent tapering cable, the 
following conditions are to be met: 

F{Z;K) = exp{K{Z-Z )) (47) 



F(Z; K) = exp{KZ) (48) 
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Then [14, equation 4] gives 



n 



3/2 



1 + 



dr 
dx 



1/4 



= (r f 2 n F(Z; K) 



(49) 



where no is the number of branches at x = 0; F(Z; K) is the 
geometric ratio imposing a taper on the equivalent cable and 
K < is the rate of taper. 



dA 

dx 



oc F(Z;K) 



dZ 

dx 



(50) 



If ^ = and every branch at any given x or Z has a different 
diameter, then for F(Z; K) = exp(KZ) [22] Then [14, equation 
6] gives : 



n(x) 



.7 = 1 L J=1 



n 



(51) 



dj is the diameter of the jth branch at distance x from soma. 
Z = when x = 0. Alternatively, if x represents actual dis- 
tance measured along successive branch points and branching 
occurs at distances = xq < ... < x v with rii branches between 
Xi and Xi + i where X{ < x < o^+i, then equation 51 becomes 
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[14, equation 7] 



.7=1 



n 



7 3/2 



3=1 



i = 0,1,. ..p 



(52) 



where (^) 2 < 1 and all branches are equal in diameter [8 p- 
156]. However, if F(Z; K) ^ 1 
then [14 equation 9] gives : 



Dtaper — D 



F(Z;K) 



2/3 



■j=i 



3/2 
ij 



2/3 



,* = 0,l,...p (53) 



From [14, equation 10] we arrive at 



A 



taper 



4i? a 



D 



taper 



1/2 



(54) 



if there is profuse branching F(Z; K) > 1 ,if there is paucity of 
branching F{Z\ K) < 1 

Jack et al[8] discusses the conditions underlying reducing a branched, 
tapering dendrite to an equivalent cylinder by taking three cases 
of taper and calculating the types of branching and values of Z 
that are possible with all three. To summarize here the types of 
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taper considered are [8, equations 7.49,7.50,7.51] 



r = ro(l — ax) (55) 



r = roexp(-ax) (56) 



r = ro(l + ax) 1 (57) 

where rg is the initial radius and a is the factor controlling the 
rate of taper with distance. If the dendrite is to be reduced to an 
equivalent cylinder the forms of branching need to be [8, equa- 
tions 7.52,7.53,7.54]: 

3 

n = uq(1 — ax)~z (58) 



n = n e 2ax ( ff 0+ } )* (59) 



/-, x5. aVn + 1 x i 

n = n (l + ax)5(— -= f- — )s (60) 

v 7 Vr 2 , + (1 + ax) 4 v 7 
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where uq is the number of branches at x = 0. 
The relationship between Z and x is an exact solution for ta- 
per described by equation (55). For the tapering conditions in 
equations (56) and (57) give approximate solutions [8, equations 
7.55,7.57,7.58]: 



Z = i((l + a 2 r$)*-(l - (1 - ax)*) (61) 
Xq a 



where 



Ao = (%?¥ (62) 



aX 



2 (e^)-l) (63) 



o 



Z« — ?— ((l + ax)i - 1) (64) 

ottAn 



3.2.3 Difference between \ocfor full morphology and equivalent cylinder: 

The equivalent cylinder formulae are strictly valid only for neu- 
rons that can be approximated as cylinders with uniform mem- 
brane resistivity. If the neuron under consideration has a den- 
dritic taper or a soma shunt, the formulas can lead to errors 



23 



as there can be violations of the equivalent cylinder approxima- 
tions. Even if these assumptions are satisifed, the accuracy of 
the parameters used in the formulae are also crucial [5] . Schier- 
wagen[27] discusses the many situations in which neurons do 
not obey all the assumptions required before reduction to an 
equivalent cylinder. Spinal motoneurons and superior colliculus 

3 

output neurons obey the rule but many others do not. One of 
the alternate models suggested is the morphology based branch- 
ing cable model. It only expects that the condition requiring the 
same boundary condition at all branches be adhered to. 
Given that most neurons do not follow the conditions in the 
equivalent cylinder approximations, the Xdc values calculated 
by the equivalent cylinder approach can at best be regarded as 
approximations. On the other hand, calculating the Xdc for the 
entire morphology of a complex neuron can be computationally 
expensive. Which approach is used may depend on weighing the 
errors versus efficiency needed in the particular problem. 

3.3 Xac Nontapering: 

This problem can be solved without the use of equivalent cylin- 
ders and with the use of equivalent cylinders. 
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3.3.1 Without equivalent cylinders : 

Calculate Xac branched for each branch using equation 28 vary- 
ing d, Ra and Cm where necessary. 

3.3.2 With equivalent cylinders : 

Calculate D from equation 43. Then substitute that in equation 



3.4 Xac tapering: 

This problem can be solved without the use of equivalent cylin- 
ders and with the use of equivalent cylinders. 

3.4.1 Without equivalent cylinders : 

Calculate XACtapered branched for each branch using equations 
37 or 40. 

3.4.2 With equivalent cylinders : 

Calculate Diaper from equation 53 and substitute it in equa- 
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(65) 



25 



tions 37 or 40. 



4 Space constant - chemical: 

Calcium concentration in dendrites and spines is a function of 
diffusion, buffering and pumping. Under certain conditions this 
can be reduced to become similar to the linear one dimensional 
cable equation [30]. Thus quantities like space constant, time 
constant and input resistance can be defined for the reaction 
diffusion equation too. We are looking at the diffusion of calcium 
ions in a cylinder after the influx of the calcium current Ica(x, t) 
across the membrane. The radial components of diffusion are 
neglected so we get one dimensional flow. The calcium ions while 
diffusing bind to various buffers and can be pumped out of the 
cylinder. The buffer itself can diffuse with a diffusion coefficient 
Db- The following equations give the change in concentration 
of calcium and bound calcium buffer. The one- dimensional 
diffusion equation is : 



dC(x, t) n d 2 C(x,t) 2 

ST = D ^r~ + r l(x ^ (66) 
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where C(x,t) is the concentration of Ca 2+ in \iM at time t 
and position x in response to the applied current density i(x,t) 
fA/fim 2 ,D is the diffusion constant in fim 2 /msec and the radius 
of the cable is r /im. A diffusable buffer is incorporated which 
shows second order kinetics : 



Ca 2+ + B^M (67) 



where Ca 2+ is the free calcium in fiM, B is the free buffer in 
\±M and M is the bound buffer in jiM. f in msec -1 and b in 
fiMmsec^ 1 are the rate constants. Incorporating this into the 
diffusion equation we can get the following : 



dM d 2 M 

-to-Ih-^ + fCB-m (69) 



B T = M(x,t) + B{x,t) (70) 

where is the diffusion constant of both the free and the bound 
buffer in fim 2 /msec and Bt is the concentration of the total 
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buffer. Under certain limiting conditions of low [Ca 2+ ] and fast 
kinetics of buffering compared to kinetics of diffusion, the above 
sets of nonlinear equations can be reduced to a single linear 
equation similar to the cable equation [30] . When an externally 
applied point source current term is used I(x,t), the resulting 
equation is [30, equation 13] : 

= r(^^)g - P m C(x,t) + K Mx ,t) 

(71) 

where (3 = P m is the membrane pump parameter at fim/msec 
and Koq is the constant of proportionality. This can be compared 
to the cable equation : 



av( x ,t) _ j_a?v _ i 
Cm at -2R i dx^ R m V(x ' t) + r,J (x ^> (72) 



where V(x, t) is the voltage, C m is the capacitance, Ri is the 
axial resistance, R m is the membrane resistance, Roo is the input 
resistance. Both equations are similar as can be seen by [30, 
equation 15] : 

R- 1 P m , C m r^±A Rr 1 + /3D b ; R^ (73) 
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Under this condition the space constant can be defined as [30, equa- 
tion 16]: 



Iannela and Tanaka[6] discuss in detail the general form of a 
nonlinear cable equation with calcium which is then reduced to 
a linear cable equation. The sodium calcium exchanger is ex- 
cluded and the voltage dependent calcium channels are assumed 
to be independent of calcium. This along with the rapid buffer 
approximation where [Ca]i <C gives rise to a pump which 
extrudes endogeneous calcium in a linear fashion. The buffer is 
also assumed to reach equilibrium much faster than the diffu- 
sion of calcium. This leads to a closed form relation between 
the rates of change of [Ca]i] and M. Thus the reaction diffusion 
equation can be reduced to the chemical cable equation (63). 

5 Discussion 

In this paper an attempt has been made to describe lambda ( 
also called space constant/ variable space constant/ length pa- 
rameter) under different conditions. Intuitively, the notion of a 
space constant or a variable space constant is more obvious in 




(74) 
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the case of a passive dendrite and often the equation to esti- 
mate this for passive cases is used for active cases too [26]. As 
described earlier, the space constant in a cable ( infinite, semi 
infinite and finite) gives an indication of the extent of voltage 
decay along the cable. The larger the diameter, the greater the 
passive spread of voltage. It also influences the summation of 
synaptic inputs. Spatially separated synaptic inputs in smaller 
diameter cables will sum differently than those in larger diameter 
cables. Once again it is important to distinguish between mor- 
phological length and the electrical length or electrotonic length 
of the dendrite. Thus in a given neuron, there could be different 
tree sizes but they could be of similar electrotonic lengths due 
to the differences in the diameter and the resulting space con- 
stant. As discussed by [28,(Chapter-5)]," The depolarization at 
the soma in response to a given input current occuring at any 
point on any dendritic tree, depends, for a given neuron, only on 
the electrotonic distance between input and soma. The response 
is independent of which dendritic tree receives the input and is 
independent of the geometrical details of the tree that receives 
the input and the geometrical details of all the other dendritic 
trees. The magnitude of the response is inversely proportional 
to the sum of the conductances per characteristic length(space 
constant) of the dendritic trunks." 

In the case of active conduction, the greater the frequency, the 
smaller is the A and thus greater is the attenuation of volt- 
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age. Here A is a function of V and this introduces nonlinearity 
into the system. However, by making certain approximations as 
shown by [15,17,6,7] the nonlinearities can be reduced to linear- 
ities by using the ionic cable theory. This involves postulating 
the nonlinear component at discrete locations and mathemati- 
cally introducing a Dirac delta function in the cable equation. 
With this transformation, length parameter used for the passive 
case can still be used in the active case. 

Reiterating the relationship between A and input resistance it 
can be said that A is directly proportional to input resistance 
in a finite and semi - infinite cable. It is proportional to \fZ ( 
impedance) in an active cell or where current injection is sinu- 
soidal. 

Goldstein and Rail [3] have shown that diameter, tapering and 
branching are important factors in determining action poten- 
tial propagation. They defined a geometric ratio (GR) which is 

3 2 

Sd ? 2 / da where d a is the diameter of the parent branch and dj 
is the diameter of the daughter branch. For uniform cables, if 
GR = 1, there is an impedance match and propagation is not af- 
fected. If GR < 1, a favorable impedance mismatch occurs and 
action potentials propagate effectively. However if GR > 1, the 
impedance mismatch is unfavorable and action potential prop- 
agation is inefficient. Vetter et al[29] showed in a simulation 
study which isolated morphology as the only variable, that back- 
propagation of action potentials was correlated with membrane 
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area distribution in the dendritic tree and the GR at individ- 
ual branch points. The same study showed that in cells with 
elaborate dendritic trees like Purkinje cells, backpropagation is 
insensitive to the sodium channel density over the physiological 
range. However in dopamine neurons, even low sodium channel 
density leads to efficient backpropagation. They also demon- 
strate that dendritic geometry places a limit on modulation of 
backpropagation by channel density and neurotransmitters. 
Bernander et al[l] and Rapp et al[25], show in two set of simula- 
tion studies that any individual cell is dependent on the network 
activity in which it is embedded. The former study simulated a 
layer 5 cortical pyramidal cell receiving inputs from 4000 excita- 
tory and 1000 inhibitory cells firing spontaneously between — 7 
Hz. Here r m and Ri n change by a factor of 10 ( 80 — 7msec and 
110 — 14 Mohms) while the electrotonic length of the cell changes 
by a factor of 3. In the [25] studywhich modeled a Purkinje cell 
with parallel fibres, it was also seen that even at a low firing 
rate of a few Hz, the parallel fibre activity changed the mem- 
brane conductance of the Purkinje cell. The time constant r m 
and the input resistance Ri n decrease several fold while electro- 
tonic length L and the voltage attenuation factor increase signif- 
icantly. This in turn affects the spatial and temporal processing 
of individual neurons. Thus results from a slice preparation and 
that from in vivo recording can give us different answers about 
the functions of an individual neuron. In this light, the electro- 
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tonic properties of the cell are not static but dynamic, evolving 
with the background environment of the cell. 
In summary, it could be said that random synaptic inputs on 
a dendrite will cause a depolarization at the soma depending 
on the electrotonic distance in the passive case. This in turn 
is dependent on the space constant of the dendrite. In the ac- 
tive case, the propagation of action potentials depends on the 
membrane area ratio and GR in the dendritic tree. This in turn 
is dependent on the ratio of diameters of parent and daughter 
branches. In other words, geometry does play a role here. In 
the active case there is an inverse relation between A w and the 
frequency of the signal. Finally, it is important to take into ac- 
count the role played by random, background synaptic activity 
on the neuron. As the synaptic activity is enhanced ( either by 
increase in numbers of synapses, or frequency), the electrotonic 
length and effective membrane time constant are both increased 
aiding in synaptic integration. 
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(a) dendrite 




Figure 1: Comparison of cable with and without boundary condition 
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